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Abstract. Non-vacuum static spherically-symmetric solutions in Palatini f(R) 
gravity are examined. It is shown that for generic choices of f(R), there are 
commonly-used equations of state for which no satisfactory physical solution of 
the field equations can be found within this framework, apart from in the special 
case of General Relativity, casting doubt on whether Palatini f(R) gravity can be 
considered as giving viable alternatives to General Relativity. 



The search for theories of gravity which can serve as an alternative to General 
Relativity (GR) has received a powerful stimulus from current developments in 
observational cosmology. If one accepts GR as the correct theory of gravity, then it 
would seem that the energy density of the universe must be dominated at the present 
time by a cosmological constant, or by an unknown form of energy (dark energy) 
that mimics the behaviour of a cosmological constant pQ. The various problems 
associated with this 2 have prompted many authors to question whether GR is indeed 
a completely correct theory of gravity on the classical level and to investigate possible 
alternatives which would not require the inclusion of dark energy. 

One of these alternatives is f(R) gravity in the Palatini formalism, or simply 
Palatini f(R) gravity [3], and we will be examining some of its properties in this paper. 
We first recall how this class of theories comes about as a generalization of GR. The 
Einstein equations can be derived from the Einstein-Hilbert action not only by the 
standard metric variation but also by making independent variations with respect to 
the metric and the connections (see, for example, [1]). This is called Palatini variation 
and, if one proceeds in this way, the form of the connections is a derived property 
rather than being specified separately (e.g. by specifying the Levi-Civita connection, 
as in standard GR). In Palatini f(R) gravity, this extended form of variation is applied 
to an action which is a general function of the scalar curvature f{R). 

S= lLJ ^ x ^9^ R ) + S m{9^,^), (1) 

where R = g^R^u, g is the determinant of the metric g^, Sm is the matter action 
and ip collectively denotes the matter fields. (We are using, throughout, units in 
which c = G = 1.) Note that here R^ v is not the Ricci tensor of the metric g^, but 
is constructed from the independent connections r „„. 
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Independent variation with respect to the metric and the connections gives 

F(R)R^-^f(R)g^ = 8TrT^, (2) 

V a W=gF{R)<T) = 0, (3) 



where F(R) = df/dR, T^ v = —2(—g)~ 1 / 2 5SM/Sg fJ ' 1 ' is the usual stress-energy tensor 
of the matter and V M is the covariant derivative with respect to the connections r A Ml/ . 
Setting f(R) — R leads to the standard GR Einstein equations while © in this case 
incorporates the definition of the Levi-Civita connection. Note that in order to derive 
([2]) and ([3]), one assumes that the matter action does not depend on the independent 
connections [see 0]. Using Palatini variation while allowing the matter to couple to 
the independent connections leads to metric-affine f(R) gravity 5 which is a different 
theory with enriched phenomenology [H [7] . 

Since it was shown that Palatini f(R) gravity might be able to address dark- 
energy problems [5] , a number of studies have been made of its cosmological aspects 
and of its consistency with cosmological constraints [9] as well as of the Newtonian and 
post-Newtonian limits and consistency with Solar System constraints [10] . However, 
the question of obtaining solutions describing stars and compact objects has received 
only a small amount of attention so far. In this paper, we focus on the problem of 
finding consistent solutions for static spherically-symmetric stars when f(R) =/= R. 

We first note that doing this is helped by the fact that Palatini f(R) gravity 
retains a useful characteristic of GR: the exterior spherically symmetric solution is 
unique (Birkhoff 's theorem). To see this one must take the trace of ([2]l. 

F(R)R-2f(R) = 8ttT, (4) 

where T = g^T^. For a given f(R), this is an algebraic equation in R and therefore 
it can be solved to give R as a function of T. We will not consider cases where this 
equation has no roots, since those do not give viable classical gravity theories 
Equation (j4|) also implies that if T = 0, R must be constant. If we denote the value 
of R when T = by Rq and insert it into ([3]) , this equation reduces to the covariant 
conservation of g^ u , fixing the connection to be the Levi-Civita one. This will be the 
situation in vacuum and in this case {2J reduces to 

R„ v - A(JZo) fl/u , = 0, (5) 

where R^ is now indeed the Ricci tensor of the metric and A(Rq) = i?o/4- According 
to whether Ro is zero or not, which of course depends on the choice of f(R), the 
theory reduces in vacuum to GR without or with a cosmological constant. The vacuum 
spherically symmetric solution is unique in either case, being either Schwarzschild or 
Schwarzschild-(anti-)de Sitter. 

Having determined an exterior solution, we then need to find an interior solution 
and perform a matching between the two. Recently, the generalisation of the 
Tolman-Oppenheimer-Volkoff (TOV) hydrostatic equilibrium equation for Palatini 
f(R) gravity has been derived [12] . opening the way for finding equilibrium interior 
solutions. We will consider here the matching of such interior solutions with exterior 
ones. 

We begin by reviewing the formulae that we will need. Denoting the Ricci scalar 
of g^v by R = g^R-nv and setting G M „ = — g^R/2, the field equations (|2|) and ([3]) 
can be rewritten as a single one 

= y t ^~ ( R -j)+j (v,*v„- g»vU) F - 
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^((V,,Fi(V,F)--,,,,,(VFr 



(6) 



where is the covariant derivative with respect to the Levi-Civita connection of 
and □ = t/^V^V,,. To arrive at this equation, one has to solve ([3]) for r*,, insert 
this into {2J and express the resulting equation in terms only of metric quantities (for 
an alternative method, see [6]). 

Using the static spherically symmetric ansatz 



ds 2 = -e A(r Ut 2 



e B(r) dr 2 



2 dfr 



(7) 



in ([5]), considering perfect-fluid matter with T M „ = (p + p^u^Uy + pg flt s (where p is the 
energy density, p is the pressure and u M is the fluid 4- velocity) and representing d/dr 
with a prime, one arrives at the equations 



A' 



B' = 



8nGrp 



— 8nGrp 
F ' 



Defining m to t{r) = r(l 
TOV equations [12]: 

p' = 




=r z V-^ 




(8) 
(9) 
(10) 

(11) 



)/2 and using Euler's equation, one gets the generalized 



(P + P) 



1 + 7 r(r - 2m to t) 
1 f4nr 2 p a+fl 



m to t 



m to t 



1+7 V F 



4atr p a 
F 2 



(r - 2m tot ) , (12) 



(13) 



We consider here matter which can be described by a one-parameter equation of state 
(EOS) p = p(p). When this is specified, one can in principle solve the above equations 
and derive an interior solution. However, this is hard to do in practice because the 
equations are implicit, their right-hand sides effectively including through F' and 
F" both first and second derivatives of the pressure, e.g., F' = d/dr [F(R(T))] = 
(dF/dR) (dR/dT) (dT/dp)p' . We therefore first put them in an explicit form, which 
allows us not only to solve them numerically, but also to study their behaviour at the 
stellar surface. 

Multiplying (|12p by dF j dp and using the definitions of a and 7, we get a quadratic 
equation in F' , whose solution is 



F' = 



-4rF(C -F)(r- 2m to t) + £>V2A 



where D 



r 2 (ZC-AF){r- 
±1 and where we have defined 



2m to t) 



dF , . dF dp , 
dp dp dp 

A = Fr 2 (r - 2m tot ) [8F(C - F) 2 (r - 2m tot ) - 

- C(4F - 3C) ((167Tp -FR + f)r 3 + 4Fm to t)] • 



(14) 

(15) 
(16) 
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We will now focus on polytropic EOSs given by p — npo T , where po is the rest- 
mass density and k and T are constants, noting that this can be rewritten as 
p = {p/n) 1 ^ + p/(T — 1), giving a direct link between p and p. In (|15p . we have 
written C in terms of dF/dp because this is finite at the stellar surface (r = r ou t where 
p = p = 0). In fact, dF/dp = (dF/dR) (dR/dT) {Up I dp - 1), where dF/dR and 
dR/dT are in general finite even when T = 3p — p goes to zero [check for instance the 
R 2 or 1/R models] and dp/ dp — > for p — > 0. Note also that while dp/ dp diverges 
when p — > 0, the product (p + p) dp/ dp goes to zero for p — > if T < 2. Therefore, for 
a polytrope with r < 2, C = at the surface. 

We now match the interior solution to the exterior one. For the latter, the general 
solution to ([5]) is exp(— B(r)) = £exp(A(r)) = 1 — 2m/r — R a r 2 /12, where I and m 
are integration constants to be fixed by requiring continuity of the metric coefficients 
across the surface and Rq is again the vacuum value of R. Using the definition of 
m to t(r) this gives, in the exterior, m to t(»') = m + r 3 i? /24 . Besides continuity of 
the metric, one has to impose some junction conditions for A'. The exterior solution 
evaluated at the surface gives 

2 (r* ont Ro - 12m) 
1 ° Ut ' ~ r out (Rot** - 12r out + 24m) ' [U) 

whereas the value of A'(r out ) for the interior solution can be calculated with (jSJ). 
For this we need F'(r ou t). Evaluating (|14|) at the surface, where C = p = and 
R, F and / take their constant vacuum values i?o, Fq and /o = FqRq/2, we get 
^"(r-out) = -(1 + D)F /r out , where D = Dsign(r out - 2m to t) S- Choosing 5 = 1 
implies 7 = — 1 at the surface [see (fTTj) ] giving A' — > 00 for r — > r~ ut [see (J5J], 
whereas A' keeps finite for r — > r+ ut [see (fT7)> ] . Because G M1/ involves A", this infinite 
discontinuity leads to the presence of Dirac deltas in the field equations. These Dirac 
deltas cannot be cancelled by the derivatives of F' on the right-hand side, because 
the discontinuity of F' is only a finite one, and one should therefore invoke an infinite 
surface density at r = r out . Since this is unreasonable, we focus only on D = — 1, 
for which F'(r out ) — when r — > r~ ut , making both F' and A' continuous across the 
surface. 

In order to study the behaviour of m to t at the surface, we need first to derive an 
explicit expression for F" . If we take the derivative of (fT4|) . F" appears on the left- 
hand side and also on the right-hand side [through m' tot , calculated from (Ti~3|) and the 
definition of j3, (fTTj) ]. giving a linear equation in F". The solution to this, evaluated 
at the surface, is 

F "(r ) - ^° r ° ut ~ 8mtot ^ c ' nsi 

r (rout) - -7. / x r l-LoJ 

8r out {r ut - 2m to t ) 

Evaluating a, (3 and 7 at the surface using F' = and F" given by (JTS]), and inserting 
into (Ti"3)) gives 

/ , , _ 2Fo^or-g ut + (rg ut gp - 8m tot ) g 

m totVoutJ — I6F0 ' 

For 1 < r < 3/2, C = dC/dpp 1 oc dC/dp (p + p) — * at the surface so that expression 
(|19p is finite and it even gives continuity of m£ ot across the surface. However, for 
3/2 < T < 2, C — > 00 as the surface is approached, provided that dF/dR(Ro) 7^ 

I Unlike in GR, one cannot prove that r ou t > 2mtot from (1 1 2 I t because p' is not necessarily positive, 
although one may expect r ou t > 2mtot in sensible solutions. 
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and dR/dT(T ) ^ (note that these conditions are satisfied by generic forms of f(R), 
i.e. whenever an R 2 term or a term inversely proportional to R is present). While 
m tot keeps finite [as can be shown using p ~ (r out — r) r / T- 1 ) , which can be derived by 
integrating (|12p near the surface] , the divergence of m' tot drives to infinity the Riemann 
tensor of the metric, R ilua \, and curvature invariants, such as R or R^ va ^ R^ a \, as 
can easily be checked [sj This singular behaviour would cause unphysical phenomena, 
such as infinite tidal forces which would destroy anything present at the surface [cf. 
the geodesic deviation equation]. We can then conclude that no physically relevant 
solution exists for any polytropic EOS with 3/2 < T < 2. Of course, polytropes give 
only simplified models for stars and the EOS in the outer layers is critical for the 
behaviour of m' tot at the surface in the non-GR case. One would like to use a more 
accurate EOS, but while this can give regular solutions in many cases {e.g. if p oc po 
near the surface), the existence of counter-examples is worrying for the viability of the 
theory. 

Setting aside the surface singularity, we next turn to the behaviour in the interior, 
focusing on models of neutron stars constructed using an analytical approximation to 
the FPS EOS [14 . Adding positive powers of R to the Einstein-Hilbert action produces 
significant effects for compact stars while adding negative ones predominantly affects 
more diffuse stars. Generically, though, one would expect terms of both types to be 
added if there is a deviation away from GR. Since the 1/R term commonly used in 
cosmology would have a negligible effect in the interior of a neutron star, we used here 
f(R) = R + eR 2 ■ As can be seen from ([6]), the metric will be sensitive to derivatives 
of the matter fields, since R is a function of t[]]. 

This can be seen in Figure [T] mtot, which in GR has a smooth profile, now 
develops peculiar features when dp /dp and d 2 p/dp 2 change rapidly in going from the 
core to the inner crust and from the inner crust to the outer crust. If mtot were plotted 
against the radius, these features would look much more abrupt, because they occur 
in a small range of radii close to the surface. While m to t does not represent a real 
mass in the interior, such a strong dependence of the metric on the derivatives of the 
matter fields is not very plausible and could have dramatic consequences. 

We have therefore found two unappealing characteristics of Palatini f(R) gravity 
as applied to stellar models, each of which arises because of the dependence of the 
metric on higher order derivatives of the matter field. First: whether or not a regular 
solution can be found depends crucially on the microphysics, through the EOS, with 
polytropic EOSs having 3/2 < T < 2 being ruled out for generic f(R). Second: even 
if an EOS does allow a regular solution at the surface, the interior metric depends 
on the first and second derivatives of the density with respect to the pressure, giving 
a problematic behaviour. While polytropic EOSs are highly idealized, we note that 
r = 5/3, corresponding to an isentropic monatomic gas or a degenerate non-relativistic 
particle gas, falls within the range not giving a regular solution. The fact that the 
gravity theory cannot provide a consistent description for these cases, casts doubt on 
whether it should be considered as a viable alternative to GR. 

Since the problems discussed here arise due to the dependence of the metric 
on higher order derivatives of the matter fields, one can expect that they will also 

§ This seems to have been missed in |15| . 

|| The unusual behaviour of this class of theories has been mentioned in a different context in |16| . 
However, we disagree with the claims made there about the violation of the equivalence principle, 
because they seem to be based on an ill-posed identification of the metric whose geodesies should 
coincide with free-fall trajectories. 
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" 1 1 11 1012 1013 1014 1015 

P (g/cm 3 ) 

Figure 1. Profiles of mtot (in Mq) and other associated quantities plotted 
against density in the interior of a neutron-star with central density 10 15 g/cm 3 
and p' = in the centre as required by local flatness. We have used the FPS 
EOS 1 141 and f(R) = R + eR 2 . The dot-dashed line shows mtot as calculated 
with e = 0.1 and the dotted line shows the equivalent curve in GR (e = 0); the solid 
line shows dmtot/dr (in Mq /'km) for e = 0.1 (this value is orders of magnitude 
lower than Solar System constraints 13 ). Note the bumps in the dmtot/dr 
curve resulting from rapid composition changes in the EOS (the corresponding 
features in the m to t curve for e = 0.1 are less apparent but a noticeable dip 
is seen at p ~ 10 14 g/cm 3 ). To make evident the influence of composition 
changes, we also show comparisons between the FPS EOS and a polytrope (with 
r = 4/3 and K = 10 15 cgs): the long-dashed and short-dashed curves show 
0.1 X (dp/dp) FPS /(o!p/c(p) poly t r0 p and 0.01 X (d 2 p/dp 2 ) F -p S / (d 2 p/dp 2 ) polytIopc , 
respectively. 



appear in other gravity theories having these characteristics. Any theory having a 
representation in which the field equations include second derivatives of the metric 
and higher than first derivatives of the matter fields will face similar problems because 
having a higher differential order in the metric than in the matter field is what 
guarantees that the metric depends in a cumulative way on the matter. If this is 
not the case then the metric loses its immunity to rapid changes in matter gradients 
since it is directly related to them instead of being an integral over them. 

The same problem should be expected for any theory which includes fields other 
than the metric for describing the gravitational interaction (e.g. scalar fields) which 
are algebraically related to matter rather than dynamically coupled. In this case 
one can always solve the field equations for the extra field and insert the solution 
into the field equation for the metric, inducing a dependence of the metric on higher 
derivatives of the matter fields. An example of such a theory is a scalar-tensor theory 
with Brans-Dicke parameter ui = —3/2, which is anyway an equivalent representation 
of Palatini f(R) gravity [6]. One should mention that this problem could probably 
be addressed in Palatini f(R) gravity by adding higher order curvature invariants in 
the action [e.g. f(R, R^R^)], since this would introduce more dynamics and break 
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the non-dynamical coupling between matter and the extra gravitational degrees of 
freedom. 

In conclusion, we suggest that our results cast doubt on the viability of theories 
including higher order derivatives of the matter fields in one of their representations, 
such as Palatini f(R) gravity or lu = —3/2 scalar-tensor theory. 

Note added: Prior to publication, a paper by Kainulainen et al. appeared [17j . which 
agreed with the validity of our results but criticized our interpretation of them. A 
response to this criticism, as well as an additional analysis of the issues presented here, 
can be found in [15] . 
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